Superfluid Quenching of the Moment of Inertia in a Strongly Interacting Fermi Gas 
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We report on the observation of a quenched moment of inertia as resulting from superfluidity in 
a strongly interacting Fermi gas. Our method is based on setting the hydrodynamic gas in slow 
rotation and determining its moment of inertia by detecting the precession of a radial quadrupole ex- 
citation. The measurements distinguish between the superfluid or collisional origin of hydrodynamic 
behavior, and show the phase transition. 
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Superfluidity is a striking property of quantum fluids 
at very low temperatures. For bosonic systems, impor- 
tant examples are liquids and clusters of 4 He and atomic 
Bose-Einstein condensates. In fermionic systems, super- 
fluidity is a more intricate phenomenon as it requires 
pairing of particles. Fermionic superfluidity is known to 
occur in atomic nuclei and 3 Hc liquids and it is also at 
the heart of superconductivity, thus being of great tech- 
nological importance. Recent advances with ultracold 
Fermi gases have opened up unprecedented possibilities 
to study theproperties of strongly interacting fermionic 
superfluids [if, 0] • Early experiments on ultracold Fermi 
gases with resonant interparticle interactions compiled 
increasing evidence for superfluidity 0, 0, @, @, 0, S| until 
the phenomenon was firmly established by the observa- 
tion of vortex lattices Q . 

In this Letter, we report on the manifestation of su- 
perfluidity in a quenched moment of inertia (MOI) in 
a strongly interacting Fermi gas undergoing slow rota- 
tion. The basic idea of a quenched MOI as a signature 
of superfluidity goes back more than 50 years in nuclear 
physics, where MOIs below the classical, rigid-body value 
were attributed to superfluidity The quenching of 

the MOI was also shown in liquid 4 Hc ll| and has, more 
recently, served for the discovery of a supersolid phase 
[T3 |. Here we introduce the observation of the quenched 
MOI as a new method to study superfluidity in ultracold 
Fermi gases. 

Our basic experimental situation is illustrated in 
Fig. [TJ At a finite temperature below the critical tem- 
perature T c , the harmonically trapped cloud consists of a 
superfluid core centered in a collisionally hydrodynamic 
cloud. We assume that the trapping potential is close 
to cylindrical symmetry, but with a slight, controllable 
deformation that rotates around the corresponding axis 
with an angular velocity f^trap- The nonsupcrfluid part 
of the cloud is then subject to friction with the trap and 
follows its rotation with an angular velocity f2, which in a 
steady state ideally reaches fl = fitrap- The correspond- 
ing angular momentum can be expressed as L = Of2, 
where O denotes the MOI. The superfluid core cannot 
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Figure 1: Schematic illustration of a strongly interacting 
Fermi gas in a slowly rotating trap. The normal part in the 
outer trap region rotates with a frequency Q, which in an 
equilibrium state corresponds to the rotation frequency fitrap 
imposed by the trap. The superfluid does not carry angular 
momentum and does therefore not contribute to the MOI. 



carry angular momentum, assuming that vortex nuclc- 
ation is avoided, and does therefore not contribute to the 
MOI of the system. Thus O represents the effective MOI 
of the whole system. 

The case of a rotating system in a steady state, where 
the normal part carries the maximum possible angu- 
lar momentum, allows us to distinguish the superfluid 
quenching of the MOI from another effect studied in 
Ref. The latter originates from irrotational flow 

and may also be discussed in terms of a MOI below the 
rigid-body value, but it can occur for both the superfluid 
and the collisionally hydrodynamic normal phase. 

Our measurements rely on the possibility to deter- 
mine the total angular momentum L of a rotating hy- 
drodynamic cloud by detecting the precession of a radial 
quadrupole excitation. This method is well established 
and has been extensively used in the context of atomic 
Bose-Einstein condensates [H, 15, lj| and recently also 
with strongly interacting Fermi gases [T3|. The preces- 
sion frequency can be written as f2 pr oc = £/(2© r i g ) 
where O r i g corresponds to a moment of inertia as calcu- 
lated from the density distribution under the assumption 
that the whole cloud, including the superfluid part, would 
perform a rigid rotation. Substituting 017 for L, we ob- 
tain f2 nr , 



= e/(20 rig ) O, with e/e rig = 1 for the full 
MOI in a normal system, and 0/O r i g < 1 for a MOI that 
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is quenched because of the supcrfluid core. 

The starting point of our experiments is an optically 
trapped, strongly interacting Fermi gas consisting of an 
equal mixture of 6 Li atoms in the lowest two atomic 
states [lj||2(J. The broad 834-G Fcshbach resonance Q 
allows us to control the s-wave interaction. If not other- 
wise stated, the measurements presented here refer to the 
resonance center, where a unitarity-limitcd Fermi gas is 
realized [H, 0] . The cigar-shaped quantum gas is confined 
in a far red-detuned, single-beam optical dipole trap with 
additional axial magnetic confinement. The trap can be 
well approximated by a harmonic potential with radial 
oscillation frequencies oj x = uj y m 2n x 680 Hz and an ax- 
ial frequency of oj z = 2tt x 24 Hz. The Fermi energy of the 
noninteracting gas is given by Ep = h(3Nu} x bjyUJ z ) 1 / 3 , 
where N = 6 x 10 5 is the total atom number. The Fermi 
temperature is Tp = Ep/k = 1.3 /iK, with k denoting 
the Boltzmann constant. 

Our scheme to study the rotational properties is de- 
scribed in detail in Ref. [l7|. It is based on a rotat- 



ing elliptical deformation of the trap, characterized by a 
small ellipticity parameter [13] e' = 0.1. In contrast to 
our previous work, we use a lower rotation frequency of 
f^trap = 27r x 200 Hz « 0.3 a-^. This low value allows us 
to avoid a resonant quadrupole mode excitation, which 
is known as an efficient mechanism for vortex nucleation 
22j . The excitation and detection of the precessmg 
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radial quadrupole mode also follows the procedures de- 
scribed in Ref. 



17] 



At this point it is important to discuss the conse- 
quences of residual trap imperfections, still present when 
we attempt to realize a cylindrically symmetric optical 
potential. As we showed in previous work 17[, we can 
control the ellipticity down to a level of ~1%. More- 
over, deviations from perfect cylindrical symmetry may 
occur because of other residual effects, such as corruga- 
tions of the optical trapping potential. As a consequence, 
a certain rotational damping is unavoidable, but damping 



times can reach typically one second (17|. This has two 
main effects for our observations. First, our measure- 
ments yield precession frequencies slightly below f2 pr cc. 
This is because of a delay time of 20 ms between turn- 
ing off the rotating trap ellipticity and applying the sur- 
face mode excitation. It is introduced to make sure that 
any possible collective excitation resulting from the ro- 
tating trap has damped out when the mode precession is 
measured. Because of rotational damping during this 
delay time, the measured precession frequencies Q' piec 
are somewhat below ri pl0 c. To compensate for this ef- 
fect, we determine the corresponding damping parame- 
ter k = f2 prec /f2 prec for each set of measurements, finding 
day-to-day variations with typical values between 0.85 
and 0.9. The second effect is induced by friction with 
static (nonrotating) trap imperfections when the rotat- 
ing ellipticity is applied. This leads to equilibrium values 
for Q typically a few percent below ^trap- For this sec- 
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Figure 2: Precession parameter V versus spin-up time t ap in for 
various values of the final temperature, as characterized by the 
heating parameter AT (see text). The quenching of the MOI 
shows up in the temperature-dependent saturation behavior. 
The applied timing sequence to facilitate measurements at 
constant temperature is illustrated above the graph. In these 
sets of measurements k = 0.85. 



ond effect there is no straightforward compensation, and 
it needs to be explicitly discussed when interpreting the 
experimental results. 

Thermometry is performed after the whole experimen- 
tal sequence. We damp out the rotation by stopping the 
trap rotation and keeping the ellipticity [23[ . We convert 
the gas into a weakly interacting one by a slow magnetic 
field ramp to 1132 G, and we finally measure the temper- 
ature T [17[ . Note that the isentropic conversion tends to 
decrease the temperature such that T is always somewhat 
below the temperature T at unitarity j24[. The relative 
statistical uncertainty of the temperature measurement 
is about 5% in the relevant temperature range. 

For our further discussion, we introduce a dimension- 
less precession parameter V by normalizing our observ- 
able Oproc to its maximum possible value of f2trap/2, 

^prec O ^ 



V 
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^trap 



(i) 



Values V < 1 show the presence of at least one of the 
two effects, namely the incomplete rotation of the nor- 
mal part (fi/fitrap < 1) or the superfluid quenching of the 
MOI (8/0 r i g < 1). Our basic idea to distinguish between 
these two effects is based on the fact that 0/O r i g rep- 
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resents a temperature-dependent equilibrium property, 
whereas Q/£l trap reflects the spin- up dynamics. There- 
fore an increase of V observed at a fixed temperature 
can be attributed to the latter effect. This, however, is 
experimentally not straightforward because of the pres- 
ence of residual heating. In the rotating trap we always 
observe some heating, which under all our experimen- 
tal conditions can be well described by a constant rate 
a = 170nK/s = 0.13 Tp/s. 

In a first set of experiments, to observe the spin-up 
dynamics at a fixed temperature, we apply a special pro- 
cedure based on the timing scheme illustrated on top of 
Fig. [21 Our procedure takes advantage of the constant 
heating rate a to control the final temperature of the 
gas. We apply the trap rotation in two separate stages 
of duration iheat and i S pin- In an intermediate time in- 
terval of £damp == 200 ms we damp out the rotation that 
is induced by the first stage [25j |. The angular momen- 
tum disappears, but the heating effect remains [23| . The 
second stage spins the cloud up again and induces fur- 
ther heating. When t tot = thcat + ispin is kept constant, 
we find that the total heating by the two rotation stages 
is AT = a itot- As only the second stage leads to a fi- 
nal angular momentum, V(t sp i n ) reveals the spin-up dy- 
namics. A variation of the parameter f tot controls AT 
and thus allows for a variation of the final temperature 
T = % + AT ; the temperature offset T w 0.11 Tp is 
set by the initial cooling and some unavoidable heating 
during the experimental sequence without trap rotation. 

Our experimental results for V(t sp i n ) are shown in 
Fig. [2] for four different values of the heating parame- 
ter AT /T F in a range between 0.026 to 0.104. This 
corresponds to a temperature range between about 0.14 
and 0.21 Tp, where we expect the supcrfluid phase tran- 
sition to occur. All four curves show qualitatively the 
same behavior. Within a few 100 ms, V rises before it 
saturates. It is straightforward to interpret this time- 
dependent increase of V in terms of the spin-up dynam- 
ics. The observed behavior can be well fit by simple ex- 
ponential curves (solid lines) approaching a steady-state 
value. This fit yields a time constant of ~50ms for all 
four curves and allows us to extract the different satura- 
tion values for V '. Note that a characteristic spin- up time 
of about 50 ms is consistent with theoretical predictions 
for a collisionally hydrodynamic gas 17, 26| . 

The four curves in Fig. [2] clearly show a general tem- 
perature dependence of the saturation value of V . For 
the two highest temperatures, however, both curves ap- 
proach nearly the same saturation value of ^0.95. This 
indicates that a superfluid core is very small or absent 
and that in both cases the hydrodynamic cloud rotates 
with an angular velocity Q near the maximum possible 
value of f2trap/2. For the lower two temperatures, how- 
ever, the saturation value shows a clear temperature de- 
pendence. At AT /T F = 0.026 and 0.052, we observe V 
to saturate at 0.68 and 0.81. Normalizing these values 
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Figure 3: The precession parameter V as a function of the 
rotation time t ro t (lower scale); the upper scale shows the 
corresponding temperature T. The vertical dashed lines sep- 
arate the three different stages (I-III) discussed in the text. 
The second stage II corresponds to the case where the normal 
part is fully rotating and the MOI is quenched because of the 
presence of a superfluid core. For this measurement k = 0.90. 



to the maximum observed V at higher temperatures, we 
derive quenching factors for the MOI of 0/O r ; g = 0.73 
and 0.87, respectively. 

Let us comment on the possible influence of vortices. 
We cannot exclude their presence, as their nuclcation can 
proceed not only via a resonant quadrupole mode exci- 
tation [2ll . |22| , but also via a coupling to the thermal 
cloud [271 ]. Vortices would result in additional angular 
momentum in the rotating cloud and its collective behav- 
ior would be closer to the normal case. This would tend 
to increase V at lower temperatures, and would counter- 
act the behavior that we observe. 

In a second set of experiments, we study the super- 
fluid phase transition in a way which is experimentally 
simpler, but which requires the knowledge of the spin-up 
dynamics as obtained from the measurements presented 
before. Here the trap rotation is applied without any 
interruption and we observe the increase of V with the 
rotation time t rot ; all other parameters and procedures 
are essentially the same as in the measurements before. 
In this procedure the temperature is not constant, but 
rises according to T = T + crfrot , where the heating rate 
a = 170nK/s is the same as before and Tq = Q.085Tp 
is somewhat lower because of the less complex timing 
sequence. 

Fig. [3] shows how V increases with the rotation time 
trot; the upper scale shows the corresponding temper- 
ature T. The observed increase in V generally results 
from both factors in Eq.(TjQ), the rising f2/f2trap and the 
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rising 8/8 r i g . However, both effects contribute on differ- 
ent time scales, which allows us to distinguish between 
three different stages (I-III) , as indicated in Fig. [3] For 
trot < 500 ms (stage I) the spin-up dynamics dominates 
and leads to a fast increase in V . For intermediate times, 
500ms < t rot < 840ms (stage II), V still increases, but 
at a much slower rate. Here, for a fully rotating normal 
part, the superfhhd quenching of the MOI (6/9 r i g < 1) 
is the dominant effect. For t rot > 840ms (stage III), V 
saturates and acquires a value that is very close to the 
maximum, indicating the full rotation of a gas without 
any supcrfluid component. 

For the supcrfluid phase transition, which occurs be- 
tween stages II and III, we obtain a critical tempera- 
ture of T c /Tp = 0.195(15). The conversion of this value, 
measured in the weakly interacting regime, to the ac- 
tual temperature in the unitarity limit yields values for 
T c between 0.20 Tp [HI and 0.24 T F This range 

for Tr is consistent with previous experimental results 
[IBB, El and theoretical predictions [JED]. 

The supcrfluid quenching effect in stage II can be 
quantitatively discussed in terms of the slope s = 
A(G/O r i g )/AT near the phase transition. Our data 
in Fig. [3] show an increase of the MOI according to 
s = 4.2(5) Tp 1 . This is fully consistent with the observed 
temperature dependence of the MOI in Fig. [21 from which 
we obtain s = 3.5(4) Tp 1 . In terms of the critical temper- 
ature T c this slope corresponds to s « LOT" 1 . We can 
compare this behavior with available predictions for the 
limits of BEC and BCS theory. For a BEC, the super- 
fluid quenching effect is very small [32| . For a BCS-type 
system, Refs. [33|, |34| predict a slope of ~ 2.0T < T 1 . Our 
observed behavior thus lies in between the two extremes 
of BEC and BCS. 

In conclusion, we have demonstrated the quenching of 
the moment of inertia that occurs as a consequence of su- 
perfluidity in a strongly interacting Fermi gas. This effect 
provides us with a novel probe for the system as, in con- 
trast to other common methods such as expansion mea- 
surements and studies of collective modes, it allows us to 
distinguish between the two possible origins of hydrody- 
namic behavior, namely collisions in a normal phase and 
superfluidity. 

We thank L. Sidorenkov for discussions. Wc ac- 
knowledge support by the Austrian Science Fund (FWF) 
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part 4). 
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